\(\int \frac {1}{(c+a^2 c x^2)^3 \arctan (a x)^{5/2}} \, dx\) [1071]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (verified)
   Maple [A] (verified)
   Fricas [F(-2)]
   Sympy [F]
   Maxima [F(-2)]
   Giac [F(-1)]
   Mupad [F(-1)]

Optimal result

Integrand size = 21, antiderivative size = 125 \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=-\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {4 \sqrt {2 \pi } \operatorname {FresnelC}\left (2 \sqrt {\frac {2}{\pi }} \sqrt {\arctan (a x)}\right )}{3 a c^3}-\frac {8 \sqrt {\pi } \operatorname {FresnelC}\left (\frac {2 \sqrt {\arctan (a x)}}{\sqrt {\pi }}\right )}{3 a c^3} \]

[Out]

-2/3/a/c^3/(a^2*x^2+1)^2/arctan(a*x)^(3/2)-8/3*FresnelC(2*arctan(a*x)^(1/2)/Pi^(1/2))*Pi^(1/2)/a/c^3-4/3*Fresn
elC(2*2^(1/2)/Pi^(1/2)*arctan(a*x)^(1/2))*2^(1/2)*Pi^(1/2)/a/c^3+16/3*x/c^3/(a^2*x^2+1)^2/arctan(a*x)^(1/2)

Rubi [A] (verified)

Time = 0.21 (sec) , antiderivative size = 125, normalized size of antiderivative = 1.00, number of steps used = 14, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.381, Rules used = {5022, 5088, 5090, 4491, 3385, 3433, 5024, 3393} \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\frac {16 x}{3 c^3 \left (a^2 x^2+1\right )^2 \sqrt {\arctan (a x)}}-\frac {2}{3 a c^3 \left (a^2 x^2+1\right )^2 \arctan (a x)^{3/2}}-\frac {4 \sqrt {2 \pi } \operatorname {FresnelC}\left (2 \sqrt {\frac {2}{\pi }} \sqrt {\arctan (a x)}\right )}{3 a c^3}-\frac {8 \sqrt {\pi } \operatorname {FresnelC}\left (\frac {2 \sqrt {\arctan (a x)}}{\sqrt {\pi }}\right )}{3 a c^3} \]

[In]

Int[1/((c + a^2*c*x^2)^3*ArcTan[a*x]^(5/2)),x]

[Out]

-2/(3*a*c^3*(1 + a^2*x^2)^2*ArcTan[a*x]^(3/2)) + (16*x)/(3*c^3*(1 + a^2*x^2)^2*Sqrt[ArcTan[a*x]]) - (4*Sqrt[2*
Pi]*FresnelC[2*Sqrt[2/Pi]*Sqrt[ArcTan[a*x]]])/(3*a*c^3) - (8*Sqrt[Pi]*FresnelC[(2*Sqrt[ArcTan[a*x]])/Sqrt[Pi]]
)/(3*a*c^3)

Rule 3385

Int[sin[Pi/2 + (e_.) + (f_.)*(x_)]/Sqrt[(c_.) + (d_.)*(x_)], x_Symbol] :> Dist[2/d, Subst[Int[Cos[f*(x^2/d)],
x], x, Sqrt[c + d*x]], x] /; FreeQ[{c, d, e, f}, x] && ComplexFreeQ[f] && EqQ[d*e - c*f, 0]

Rule 3393

Int[((c_.) + (d_.)*(x_))^(m_)*sin[(e_.) + (f_.)*(x_)]^(n_), x_Symbol] :> Int[ExpandTrigReduce[(c + d*x)^m, Sin
[e + f*x]^n, x], x] /; FreeQ[{c, d, e, f, m}, x] && IGtQ[n, 1] && ( !RationalQ[m] || (GeQ[m, -1] && LtQ[m, 1])
)

Rule 3433

Int[Cos[(d_.)*((e_.) + (f_.)*(x_))^2], x_Symbol] :> Simp[(Sqrt[Pi/2]/(f*Rt[d, 2]))*FresnelC[Sqrt[2/Pi]*Rt[d, 2
]*(e + f*x)], x] /; FreeQ[{d, e, f}, x]

Rule 4491

Int[Cos[(a_.) + (b_.)*(x_)]^(p_.)*((c_.) + (d_.)*(x_))^(m_.)*Sin[(a_.) + (b_.)*(x_)]^(n_.), x_Symbol] :> Int[E
xpandTrigReduce[(c + d*x)^m, Sin[a + b*x]^n*Cos[a + b*x]^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && IGtQ[n, 0]
&& IGtQ[p, 0]

Rule 5022

Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> Simp[(d + e*x^2)^(q + 1)*
((a + b*ArcTan[c*x])^(p + 1)/(b*c*d*(p + 1))), x] - Dist[2*c*((q + 1)/(b*(p + 1))), Int[x*(d + e*x^2)^q*(a + b
*ArcTan[c*x])^(p + 1), x], x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[e, c^2*d] && LtQ[q, -1] && LtQ[p, -1]

Rule 5024

Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_.)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> Dist[d^q/c, Subst[Int[(a
 + b*x)^p/Cos[x]^(2*(q + 1)), x], x, ArcTan[c*x]], x] /; FreeQ[{a, b, c, d, e, p}, x] && EqQ[e, c^2*d] && ILtQ
[2*(q + 1), 0] && (IntegerQ[q] || GtQ[d, 0])

Rule 5088

Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_.)*(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> Simp[x^m*(d +
 e*x^2)^(q + 1)*((a + b*ArcTan[c*x])^(p + 1)/(b*c*d*(p + 1))), x] + (-Dist[c*((m + 2*q + 2)/(b*(p + 1))), Int[
x^(m + 1)*(d + e*x^2)^q*(a + b*ArcTan[c*x])^(p + 1), x], x] - Dist[m/(b*c*(p + 1)), Int[x^(m - 1)*(d + e*x^2)^
q*(a + b*ArcTan[c*x])^(p + 1), x], x]) /; FreeQ[{a, b, c, d, e, m}, x] && EqQ[e, c^2*d] && IntegerQ[m] && LtQ[
q, -1] && LtQ[p, -1] && NeQ[m + 2*q + 2, 0]

Rule 5090

Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_.)*(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> Dist[d^q/c^(m
 + 1), Subst[Int[(a + b*x)^p*(Sin[x]^m/Cos[x]^(m + 2*(q + 1))), x], x, ArcTan[c*x]], x] /; FreeQ[{a, b, c, d,
e, p}, x] && EqQ[e, c^2*d] && IGtQ[m, 0] && ILtQ[m + 2*q + 1, 0] && (IntegerQ[q] || GtQ[d, 0])

Rubi steps \begin{align*} \text {integral}& = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}-\frac {1}{3} (8 a) \int \frac {x}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{3/2}} \, dx \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {16}{3} \int \frac {1}{\left (c+a^2 c x^2\right )^3 \sqrt {\arctan (a x)}} \, dx+\left (16 a^2\right ) \int \frac {x^2}{\left (c+a^2 c x^2\right )^3 \sqrt {\arctan (a x)}} \, dx \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {16 \text {Subst}\left (\int \frac {\cos ^4(x)}{\sqrt {x}} \, dx,x,\arctan (a x)\right )}{3 a c^3}+\frac {16 \text {Subst}\left (\int \frac {\cos ^2(x) \sin ^2(x)}{\sqrt {x}} \, dx,x,\arctan (a x)\right )}{a c^3} \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {16 \text {Subst}\left (\int \left (\frac {3}{8 \sqrt {x}}+\frac {\cos (2 x)}{2 \sqrt {x}}+\frac {\cos (4 x)}{8 \sqrt {x}}\right ) \, dx,x,\arctan (a x)\right )}{3 a c^3}+\frac {16 \text {Subst}\left (\int \left (\frac {1}{8 \sqrt {x}}-\frac {\cos (4 x)}{8 \sqrt {x}}\right ) \, dx,x,\arctan (a x)\right )}{a c^3} \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {2 \text {Subst}\left (\int \frac {\cos (4 x)}{\sqrt {x}} \, dx,x,\arctan (a x)\right )}{3 a c^3}-\frac {2 \text {Subst}\left (\int \frac {\cos (4 x)}{\sqrt {x}} \, dx,x,\arctan (a x)\right )}{a c^3}-\frac {8 \text {Subst}\left (\int \frac {\cos (2 x)}{\sqrt {x}} \, dx,x,\arctan (a x)\right )}{3 a c^3} \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {4 \text {Subst}\left (\int \cos \left (4 x^2\right ) \, dx,x,\sqrt {\arctan (a x)}\right )}{3 a c^3}-\frac {4 \text {Subst}\left (\int \cos \left (4 x^2\right ) \, dx,x,\sqrt {\arctan (a x)}\right )}{a c^3}-\frac {16 \text {Subst}\left (\int \cos \left (2 x^2\right ) \, dx,x,\sqrt {\arctan (a x)}\right )}{3 a c^3} \\ & = -\frac {2}{3 a c^3 \left (1+a^2 x^2\right )^2 \arctan (a x)^{3/2}}+\frac {16 x}{3 c^3 \left (1+a^2 x^2\right )^2 \sqrt {\arctan (a x)}}-\frac {4 \sqrt {2 \pi } \operatorname {FresnelC}\left (2 \sqrt {\frac {2}{\pi }} \sqrt {\arctan (a x)}\right )}{3 a c^3}-\frac {8 \sqrt {\pi } \operatorname {FresnelC}\left (\frac {2 \sqrt {\arctan (a x)}}{\sqrt {\pi }}\right )}{3 a c^3} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 0.69 (sec) , antiderivative size = 186, normalized size of antiderivative = 1.49 \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\frac {2 \left (-\frac {1}{a \left (1+a^2 x^2\right )^2}+\frac {8 x \arctan (a x)}{\left (1+a^2 x^2\right )^2}-\frac {\sqrt {2} (-i \arctan (a x))^{3/2} \Gamma \left (\frac {1}{2},-2 i \arctan (a x)\right )}{a}+\frac {\sqrt {2} \arctan (a x)^2 \Gamma \left (\frac {1}{2},2 i \arctan (a x)\right )}{a \sqrt {i \arctan (a x)}}-\frac {(-i \arctan (a x))^{3/2} \Gamma \left (\frac {1}{2},-4 i \arctan (a x)\right )}{a}+\frac {\arctan (a x)^2 \Gamma \left (\frac {1}{2},4 i \arctan (a x)\right )}{a \sqrt {i \arctan (a x)}}\right )}{3 c^3 \arctan (a x)^{3/2}} \]

[In]

Integrate[1/((c + a^2*c*x^2)^3*ArcTan[a*x]^(5/2)),x]

[Out]

(2*(-(1/(a*(1 + a^2*x^2)^2)) + (8*x*ArcTan[a*x])/(1 + a^2*x^2)^2 - (Sqrt[2]*((-I)*ArcTan[a*x])^(3/2)*Gamma[1/2
, (-2*I)*ArcTan[a*x]])/a + (Sqrt[2]*ArcTan[a*x]^2*Gamma[1/2, (2*I)*ArcTan[a*x]])/(a*Sqrt[I*ArcTan[a*x]]) - (((
-I)*ArcTan[a*x])^(3/2)*Gamma[1/2, (-4*I)*ArcTan[a*x]])/a + (ArcTan[a*x]^2*Gamma[1/2, (4*I)*ArcTan[a*x]])/(a*Sq
rt[I*ArcTan[a*x]])))/(3*c^3*ArcTan[a*x]^(3/2))

Maple [A] (verified)

Time = 0.18 (sec) , antiderivative size = 113, normalized size of antiderivative = 0.90

method result size
default \(\frac {-16 \sqrt {2}\, \sqrt {\pi }\, \operatorname {FresnelC}\left (\frac {2 \sqrt {2}\, \sqrt {\arctan \left (a x \right )}}{\sqrt {\pi }}\right ) \arctan \left (a x \right )^{\frac {3}{2}}-32 \sqrt {\pi }\, \operatorname {FresnelC}\left (\frac {2 \sqrt {\arctan \left (a x \right )}}{\sqrt {\pi }}\right ) \arctan \left (a x \right )^{\frac {3}{2}}+16 \sin \left (2 \arctan \left (a x \right )\right ) \arctan \left (a x \right )+8 \sin \left (4 \arctan \left (a x \right )\right ) \arctan \left (a x \right )-4 \cos \left (2 \arctan \left (a x \right )\right )-\cos \left (4 \arctan \left (a x \right )\right )-3}{12 a \,c^{3} \arctan \left (a x \right )^{\frac {3}{2}}}\) \(113\)

[In]

int(1/(a^2*c*x^2+c)^3/arctan(a*x)^(5/2),x,method=_RETURNVERBOSE)

[Out]

1/12/a/c^3*(-16*2^(1/2)*Pi^(1/2)*FresnelC(2*2^(1/2)/Pi^(1/2)*arctan(a*x)^(1/2))*arctan(a*x)^(3/2)-32*Pi^(1/2)*
FresnelC(2*arctan(a*x)^(1/2)/Pi^(1/2))*arctan(a*x)^(3/2)+16*sin(2*arctan(a*x))*arctan(a*x)+8*sin(4*arctan(a*x)
)*arctan(a*x)-4*cos(2*arctan(a*x))-cos(4*arctan(a*x))-3)/arctan(a*x)^(3/2)

Fricas [F(-2)]

Exception generated. \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate(1/(a^2*c*x^2+c)^3/arctan(a*x)^(5/2),x, algorithm="fricas")

[Out]

Exception raised: TypeError >>  Error detected within library code:   integrate: implementation incomplete (co
nstant residues)

Sympy [F]

\[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\frac {\int \frac {1}{a^{6} x^{6} \operatorname {atan}^{\frac {5}{2}}{\left (a x \right )} + 3 a^{4} x^{4} \operatorname {atan}^{\frac {5}{2}}{\left (a x \right )} + 3 a^{2} x^{2} \operatorname {atan}^{\frac {5}{2}}{\left (a x \right )} + \operatorname {atan}^{\frac {5}{2}}{\left (a x \right )}}\, dx}{c^{3}} \]

[In]

integrate(1/(a**2*c*x**2+c)**3/atan(a*x)**(5/2),x)

[Out]

Integral(1/(a**6*x**6*atan(a*x)**(5/2) + 3*a**4*x**4*atan(a*x)**(5/2) + 3*a**2*x**2*atan(a*x)**(5/2) + atan(a*
x)**(5/2)), x)/c**3

Maxima [F(-2)]

Exception generated. \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\text {Exception raised: RuntimeError} \]

[In]

integrate(1/(a^2*c*x^2+c)^3/arctan(a*x)^(5/2),x, algorithm="maxima")

[Out]

Exception raised: RuntimeError >> ECL says: expt: undefined: 0 to a negative exponent.

Giac [F(-1)]

Timed out. \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\text {Timed out} \]

[In]

integrate(1/(a^2*c*x^2+c)^3/arctan(a*x)^(5/2),x, algorithm="giac")

[Out]

Timed out

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{\left (c+a^2 c x^2\right )^3 \arctan (a x)^{5/2}} \, dx=\int \frac {1}{{\mathrm {atan}\left (a\,x\right )}^{5/2}\,{\left (c\,a^2\,x^2+c\right )}^3} \,d x \]

[In]

int(1/(atan(a*x)^(5/2)*(c + a^2*c*x^2)^3),x)

[Out]

int(1/(atan(a*x)^(5/2)*(c + a^2*c*x^2)^3), x)